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S ' Abstract 

''O ■ Using the path integral approach, we provide an exphcit derivation of the equa- 



^ \ tion for the phase boundary for quantum Josephson junction arrays with offset 

fj ■ charges and non-diagonal capacitance matrix. For the model with nearest neigh- 

bor capacitance matrix and uniform offset charge q/2e — 1/2, we determine, in 
the low critical temperature expansion, the most relevant contributions to the 
►^ \ equation for the phase boundary. We explicitly construct the charge distribu- 

(y-^ ■ tions on the lattice corresponding to the lowest energies. We find a reentrant 

OO \ behavior even with a short ranged interaction. 

A merit of the path integral approach is that it allows to provide an elegant 

derivation of the Ginzburg-Landau free energy for a general model with charge 

frustration and non-diagonal capacitance matrix. The partition function factor- 

O^ \ izes as a product of a topological term, depending only on a set of integers, and 

a non-topological one, which is explicitly evaluated. 



Ch ! 1 Introduction 

O 

^ ' Josephson junction arrays (JJA) and granular superconductors, namely sys- 

tems of metallic grains embedded in an insulator, become superconducting in 
k>( \ two steps ||l[. First, at the bulk critical temperature each grain develops a su- 

'V^ i perconducting gap but the phases of the order parameter on different grains 

C^ I are uncorrelated. Then, at a lower temperature T^, the Cooper pair tunneling 

between grains gives rise to a long-range phase coherence and the system as a 
whole exhibits a phase transition to a superconducting state. 

^Present address: S.I.S.S.A., Via Beirut 4, 34814 Trieste, Italy 



The phase transition is governed by the competition between the Joseph- 
son tunneling, characterized by a Josephson coupUng energy Ej g], and the 
Coulomb interaction between Cooper pairs, described by a charging energy 
Ec 1^, |4). In classical junction arrays the Josephson coupling Ej is dominant, 
the transition separates a superconducting low temperature phase from a normal 
high temperature phase. When Ec is comparable to Ej (small grains) charging 
effects give rise to a quantum dynamics. The grain capacitance is small, so that 
the energy cost of Cooper pair tunneling may be higher then the energy gained 
by the formation of a phase-coherent state. Zero point fluctuations of the phase 
may destroy the long range superconducting order even at zero temperature 
(see for example [Q). 

Within the framework of the BCS theory, Efetov [g derived an effective 
quantum Hamiltonian in terms of the phases ipi of the superconducting order 
parameter at the grain i, and their conjugate variables ni number of Cooper 
pairs. Efetov's Hamiltonian for the quantum phase model reads 

ij <ij> 

Qi = 2eni [Vi,ni] = i (Jy , 

where Qi is the excess of charge due to Cooper pairs (charge 2e) on the site i of a 
square lattice in D-space dimension and Cy is the capacitance matrix describing 
the electrostatic coupling between Cooper pairs. The diagonal elements of the 
inverse matrix C^^ provide the charging energy: Ec = e'^C^^/2 = e^/2Co, 
where Cq is the self-capacitance. 

The superconductor-insulator transition depends crucially on the spatial di- 
mensionality D. For D — 1 there may exist also other phases 0. For D — 2 
the system exhibits a richly structured phase diagram (see for example W 0]). 
In higher dimensions it is believed that the mean field theory analysis provides 
qualitatively correct results. 

It is relevant to understand how the transition from insulator to superconduc- 
tor depends on the relevant constitutive parameters - such as the capacitances 
of, and between, the junctions - as well as on external parameters such as the 
temperature, offset charges and external magnetic fields. 

Much work has been done to study the phase diagram of quantum JJA, in the 
T/Ec-Ec/Ej plane § . The analysis uses mean field theory |, |, |ll|, |l|, ^ 
|lj, |l^ as well as the renormalization group approach [|l6[ |l^ . There is the claim 
that the phase diagram -under suitable circumstances- may exhibit a reentrant 
character with the superconducting phase existing between upper and lower 
critical temperature |^, |l^, |^ . In [g|, |l^ the influence of the Coulomb energy 
on the transition temperature was investigated for a model with a diagonal 
capacitance matrix. The effects of off-diagonal terms in the charging energy were 
investigated by several authors within the mean field theory approximation p, 
|ll|, |l2[ ^ |20[ |2|] : while it is widely believed that the nearest neighbors 



interaction enhance the transition temperature Tc by lowering the energy cost 
for a Cooper pair to tunnel from one neutral grain to another p2[, there is 



still some dispute on whether there is a reentrance or not for models with non- 



diagonal capacitance matrix |12, EOl pl[ . 

It is relevant for physical applications to consider the effect of a background 
of external charges on the superconductor-insulator transition of a quantum 
J JA iQ, ^, |2^, ^ . Such an offset of charges might arise in physical systems as 
a result of charged impurities or by application of a gate voltage between the 
array and the ground. In the former case offset charges are distributed randomly 
on the lattice while in the latter case they play the role of a chemical potential 
for charges. They might be regarded as effective charges ^i on the sites of the 
lattice. When ^i 7^ 2e the offset charges cannot be eliminated by Cooper pair 
tunneling. 

Offset charges frustrate the attempts of the system to minimize the energy of 
the charge distribution of the ground state. Consequently the charging energy 
of any excitations is smaller compared to the unfrustrated case and supercon- 
ductivity is enhanced. With offset charges the Hamiltonian (||) becomes 

ij <ij> 

In order to study the effect of charge frustration on the phase diagram of the 
system described by the Hamiltonian (||), it is our purpose to revisit the mean 
field theory of quantum JJA using the path-integral method. The approach 
uses the Hubbard-Stratonovich pj] representation for the partition function 
in terms of coarse-grained classical local variables ipi for which the effective 
action is computed p9[ | . We find a reentrant behavior for models with a nearest 
neighbor capacitance matrix and a uniform offset charge qi = e, even if the 
interaction among grains is short ranged. We find analytically the equation 
which determines the critical temperature as a function of Ej/Eq- This allows 
us to analyze the low temperature limit of the theory and to find the regimes 
in which a reentrant behavior might be observed. 

In section 2 we review the self consistent mean field theory approximation 
within the Hamiltonian formalism for quantum JJA. We study the eigenvalue 
equation of the mean field Hamiltonian with diagonal capacitance, and uni- 
form offset charge qi = e showing explicitly that at zero temperature there is 
superconductivity for all values of a = zEj/AEc- 

In section 3 we use the coarse grained approach to compute the Ginzburg- 
Landau free energy for quantum JJA with charge frustration and a general 
Coulomb interaction matrix. The path integral providing the phase correlator 
needed to investigate the critical behavior of the system, is explicitly computed. 

In section 4, from the Ginzburg-Landau free energy, we derive, within the 
mean field theory approximation, the analytical form of the critical line equa- 
tion. The phase diagram is drawn in the diagonal case for a generic external 
charge distribution. We then analyze the low temperature limit of a system 
with nearest neighbor interaction matrix and find a reentrant behavior when a 
uniform background of external charges <7i = e is considered. 



Section 5 is devoted to some concluding remarks. The appendices contain 
the derivation of some relevant formulas of the main text. 



2 Self-consistent mean field theory in the Hamil- 
tonian approach 

Mean field theory for quantum JJA with diagonal capacitance matrix was first 
used by Symanek [p|. The approximation consists in replacing the Josephson 
coupling of the phase on a given island i to its neighbors by an average coupling 
so that: 

Ej 2_. cos(iy9i — ipj) = zEj < cosiy9 > 2_. cos(^i. (3) 

<ij> i 

In (H) z is the coordination number; it is assumed also that < cos ip > does not 
depend on the island index i and the choice < sin(pi >= 0, which provides a 
real order parameter, has been made. 

Within the mean field approximation the Hamiltonian (nh becomes 

-^™/ " 2 ^ Cr^^QiQi - zEj < cos ip>^ cos ipi (4) 



and the order parameter < cos ip > is evaluated self-consistently from (|j) . For 
a diagonal capacitance matrix (Cy = C'o'^ij) niean field theory computation are 
very simple since (^ describes on each site a quantum particle in a periodic 
potential cos ipi [g| . 

The eigenfunction of the array is a product of eigenfunctions ■(/'n(v) describ- 
ing the individual islands and satisfying the Mathieu equation pq| 

< COS(p > COSip^tpnif) ^ —^tpni'P) (5) 



dip^ iEc yj-rny-r, ^^^ 

with periodic boundary conditions tpniv) = '0n(<P + 27r). 

It is well known that the Mathieu equation admits also antiperiodic solutions, 
ipni'p) — —4'ni'P + 27r) (see appendix A). If both periodic and antiperiodic 
solutions are used, the general solution of (||) does not have a definite periodicity 
and, consequently, the charges ni take continuous eigenvalues. Such continuous 
eigenvalues are expected to be relevant in the description of continuous flows 
of currents due for example to ohmic shunt resistances |26| |2^. Although the 
superposition of both periodic and antiperiodic solutions yields to a reentrant 
behavior even in the unfrustrated dissipationless diagonal model ||l^, |ri|, p9| , 
this superposition is not allowed in describing physical situations in which the 
only excitations are Cooper pairs of charge 2e |^, |l^, Q . Thus the use of both 
periodic and antiperiodic solutions does not have physical significance in the 
models considered in this paper. 

The mean field self-consistency condition gives 

< COS 99 >= -^ ^^ ^_f3E^ (6) 



with /3 = 1/kBT. For high temperatures or low Ej only the solution < cos ip > = 
exists and there is not superconductivity. For low temperatures or high Ej 
< cosip >y^ and the system as a whole behaves as a superconductor. 
From (pf) one gets the equation for the critical temperature H 
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with y = kg/TcEc and a — zEj/AEq- 
In fig.l we plot Tc as a function of a. 
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Figure 1: Phase diagram for the diagonal model without charge frustration. 



If one considers a diagonal capacitance matrix and uniform offset charges of 
magnitude e on each site (qi/2e = 1/2), the Hamiltonian reads 



Mean field theory of this model leads to a Schrodinger equation of the form 

2 



dip'^ 2e dip \ 2e 



a < cos p > cos p) 



■0«(V) 



En 

4E^ 



^»(^). (9) 



Redefining the phase of the wave function as 

il>n{p) = e'*A'^p„(<^) 



reduces to a Mathieu equation for Pn{f) 
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(10) 



with A„ = En/Ec and v = -zEj < cos(^ > /2Ec- Equations (|8|) , @ , (|lO|) lead 
to the following modification of (^ [see Appendix A] : 



(11) 



which - at variance with the unfrustrated model - exhibits superconductivity 
even for infinitesimal values of a. This feature is shown in fig. 2 which also 
shows the absence of a reentrant behavior at low T. 
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Figure 2: Phase diagram of the diagonal model with half-integer charge frustra- 
tion. 



For frustrated models with non diagonal capacitance matrix, the self-consistent 
mean-field theory approximation becomes very cumbersome and one should re- 
sort to the more powerful functional approach. A reentrance at a low Tc is 
expected at least when the interaction between grains is long-ranged p2, 0| . 



3 Ginzburg-Landau free energy 

The partition function for the frustrated off-diagonal model is given by 

n 

where H is given in (0) and the sum is extended only to states of charge 2e and 
thus with definite periodicity. 

In the functional approach Z reads 

Z = I n ^^i ^^P { - / ^tLe (^¥>i(r), ^(r)) | (13) 

where the Euclidean Lagrangian Le can be derived from 



1/ h''^ 



^ ^ ij ^ ^ i <ij> 

by replacing it/h -^ r. The path integral that one should compute is then given 
by: 

i *- ^ ij i ij 

(15) 
where —oo<ipi< -|-cx), ^i{0) = '^i{(i) + 27rni and 7ij = 1 if i,j are nearest 
neighbors and equals zero otherwise. The integers ni appearing in this boundary 
conditions take into account the 27r-periodicity of the states ^„ appearing in 

®- 

In order to derive the Ginzburg-Landau free energy for the order parameter, 

it is convenient to carry out the integration over the phase variables by means 

of the Hubbard-Stratonovich procedure [p4|: using the identity 

e-^+rj ^ de^ j ^ ^,^_ g-^+r- V-j+v-^+ J (16) 

i 

the partition function may be rewritten as 

i 

where the effective action for the auxiliary Hubbard-Stratonovich field -0;, Ssff [ip] i 
is given by 



?i^//M = -log{/n^^i expj^ dr[-i^Cyg|^ 



2e 



+^E(9i£-V'ie'^'-V^re~^''')]}}- (18) 

The Hubbard-Stratonovich field ^i may be regarded as the order parameter for 
the insulator-superconductor phase transition since it turns out to be propor- 
tional to < e*'''' > , as it can be easily seen from the classical equations of motion. 
From (|l8|) the Ginzburg-Landau free-energy may be derived by integrating out 
the phase field (jji. 

Since the phase transition is second order [g8| , close to the onset of supercon- 
ductivity, the order parameter ipi is small. One may then expand the effective 
action up to the second order in ipi, getting 

rP rP 

SesM = SEffm + dr dT'GrsiT,T')MrWs{r') + ■■■ (19) 
Jo Jo 

where G^s is the phase correlator 

s'SEffm 



Grs(T,T') = 



= (e'Vr(r)-i<Ps(r')^g_ (20) 



The partition function (O), can be written as 

Z= /"jldV-idV^re-^I^l (21) 



where F[ip] is the Ginzburg-Landau free energy; due to ( ]l8| , |l9| ), up to the second 
order in ipi, one has 

F[^] ^ f dr f dT'^^rM[7i7''5(r- r') - Gij(r, r')]^j(r'). (22) 

Jo Jo jj 

We shall now compute the phase correlator Grs by evaluating the expecta- 
tion value in (p^) by means of the path integral over the phase variables (pi (r) . 
In performing this integration one should take into account that the field con- 
figurations satisfy 

ipiiP) - ipi{0) = 2Tmi . (23) 

For this purpose it turns out very convenient to untwist the boundary conditions 
by decomposing the phase field in terms of a periodic field (/)i (t) and a term linear 
in T which takes into account the boundary conditions (p3|); namely, one sets 



Viir) = 0i(r) + -^mr , (24) 

with (f>i{f3) = (l)i{0). Summing over all the phases (piir) amounts then to inte- 
grate over the periodic field (f>i and to sum over the integers ni. As a result the 



phase correlator factorizes as the product of a topological term depending on 
the integers rii and a non-topological one; namely 

Grs(r;T') = — 



Sfn. 



e p ^ 



/i5</.iexp|/o'^dr(-iCy^^) 



(25) 



Eij"if!2Cij"i"j+X;i2i^/5|^ni ^ 

After a lengthy computation, the first (non-topological) factor appearing in the 
l.h.s. of equation (^5|) has the following simple expression [see Appendix B]: 



6r. exp|-2e2C-/ (^|r - r'| - ^^^) }■ (26) 

The sum over the integers in the topological factor in (£3) is done by means of 
the Poisson resummation formula 

I detGI^ y^e~'^^"^'^^''^'-'^"^'^^j _ y^g-7rmi(G"^)ijmj-27rimiai 
[m] [nii] 

Thus eq.(p3) becomes 



Grs(r,r')-^rse-2^^"l^-^'l • 



(27) 



Z^[ni] ^ 

with rii assuming all integer values and X^fn-i being a sum over all the configu- 
rations. 

By means of a Euclidean-time Fourier transform, the fields fAi are written as 

Mr) ^ ^Y.^ii^^^y""^ ' 

where w^ are the Matsubara frequencies. As a consequence, the phase correlator 
Gij can be expressed as 

Gij(r;r') = i^Gij(c.^;c.^Oe^"-^e-^-' . (28) 



From (27) one can show that Grs(wp; tj' ) is diagonal in the Matsubara frequen- 



cies and can be written as 

Grs(t^p; W^) = Gr(w^) • Srs ' S{uJf^ + W^/ ) (29) 



with 

GrK) = Y—- 77 ■ (30) 



In (|3l| ) Zo is given by 

[ni] 

with C/ij = C;T^, -Ec = e^CrrV2 and y = ksTc/Ec- In terms of Matsubara 
frequencies the Ginzburg-Landau free energy ( p2[ ) becomes 

^jK) • (31) 



^m - ^ E^i K) [ J-7i7' - Gi{Lu,)S,i 



This is our starting point for any analysis of the phase boundary between the 
insulating and the superconducting phases in JJA with arbitrary capacitance 
matrix and with charge frustration. 

4 Mean field theory analysis 

In the following we shall derive the equation determining the phase boundary 
in the plane {a,KBTc/Ec), in mean field theory and for a system with arbi- 
trary capacitance matrix and a uniform distribution of off-set charges. For this 
purpose it is convenient to expand the fields ipi{ujf^) and Gi{ujfj,) in terms of the 
vectors of the reciprocal lattice q. One has 

V^iK) = ^E^<i(^M)e''^-' (32) 

q 

^'(^/^) = ^Y.Gciiu^'^y''' ■ (33) 



q 



Moreover 



y,' = ^T.^.'^''-'-'- (34) 

q 

where 7"^ is the inverse of the Fourier transform of the Josephson coupling 
strength 7ij which equals 1 for i,j nearest neighbors and otherwise. As a 
consequence 

-i_ 1 



7q 



^. 



g-iqp 



where p is a vector connecting two nearest neighbors sites. Expanding in q one 
gets 



10 



where a is the lattice spacing and z the coordination number. The first term in 
(pq) provides the mean field theory approximation which, as expected, is exact 
in the limit of large coordination number. 

The Ginzburg-Landau free energy (pl|), reads 

/xqq' 

Using (pq) and keeping only terms of zero-th order in w^ and q one obtains the 
mean field theory approximation to the coefficient of the quadratic term of F 



—y 

[iN ^ 



^-^°(«) 



l^qK)!' • (36) 



The equation for the phase boundary line then reads as 

l = z^Go(0) (37) 

with 

^oW^^E^rW. (38) 

r 

Equation (|37| ) determines the relation between Tc and a at the phase boundary. 
For a uniform distribution of offset charges eq. ( p7| ) simplifies further since in 
Gr does not depend on r. As a consequence, the phase boundary equation 
becomes 

^ e-^^>i^<"'+*''"^ + *^ 1 

l = a-Y j-r (39) 

U l-4Ej^(ni + e^ Z, 

with 

["i] 

In the following we shall derive the physical implications of (^9|) in a variety 
of models describing JJA. 

4.1 Self-charging model 

For a diagonal capacitance matrix, U-^ = 5i^Uo, one singles out only the self- 
interaction of plaquettes. This case was already analyzed in section 2 within the 
approach of self-consistent mean field theory. As a check of the path integral 
approach we shall show that one is able to reproduce the same results from 
eq.®. 

In the diagonal case eq.(p9|) becomes 



n 



(n+ 



9 ^2 



y e-^(™+*) 
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i-4("+^)^ 1 (40) 
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Figure 3: Phase diagram for several values of /j. 



Since n is an integer ( [fO| ) is invariant under the shift ^ -^ ^ + 1- ^'^^ q = 
eq. (|4Cl| ) reduces to (0) . From figure 1 one readily sees that there is no supercon- 
ductivity for a < 1. Due to the periodicity of (|o|) this holds for any integer 
q. For q/2e equal to 1/2 one gets equation (|ll|). From figure 2 one sees that 
superconductivity is attained for all the values of a, since the superconducting 
order parameter at zero temperature is different from zero. 

For the self-charging model the system exhibits superconductivity for all the 
values of a also if the distribution of offset charges is such that integer and 
half-integer charges coexist on the lattice. If one denotes by /o the fraction of 
integer charges and by /i = 1 — /o the fraction of half-integer charges, eq.(^0|) 
implies that 



a = ./o 



Er 



E. 



l-4(n+i)2 



Er 



Er 



t(m+i)2 



In fig. 3 we plot Tc as a function of a for several values of /g. As expected 
superconductivity is enhanced as /i increases. 

4.2 Models with non-diagonal capacitance matrix 

In pO| Fishman and Stroud, using a low temperature expansion, determined 
Tc as a function of a for models with non diagonal interaction matrix without 
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considering the effect of offset cliarges. They did not find signs of normal state 
reentrance for nearest neighbor interaction matrix models in which only the 
diagonal interaction matrix element Uqo and the nearest neighbor interaction 
matrix element C/op = 9Uoo are nonzero. This can be seen from the expansion 
of the critical line eq.(p9[) for g = and small critical temperatures: 

Reentrant behavior is possible (l2| for 9 > 9c — ,.j, „ when the coefficient of the 

exponential e~^'^ is negative; in fact, the phase boundary line a — a{Tc) first 
bends to the left due to the negative coefficient of e^'^^^ and finally, when the 
critical temperature is high enough, bends to the right, favoring the insulating 
phase. 

As evidenced by Fishman and Stroud |2C|, the regime of physical interest 
is 9 < -; namely, when the capacitance matrix is invertible. Reentrance is 
possible only in one dimension {9c = l/vTO < 1/z — 1/2); in higher dimensions 
reentrance occurs only when the electrostatic interaction is long ranged EQI . 

If there are half-integer offset charges on the sites of a square lattice, our 
analysis shows that the equation for the critical line is 

^^Z^M (41) 

In eq.(pl|) appears the expression 

ij 

which is the electrostatic energy of a generic charge distribution on the lattice. 
Denoting with n? and nl the charge distributions of the two lowest lying en- 
ergy states and with E'^ and E^ the corresponding energies, the low temperature 
expansion of eq. (^ yields 



(43) 



Independently on the explicit form of Uij, E[ni\ reaches its minimum value 
when (nP-^i) = ±i(-l)*i+i2-H- +«d .^^jth ij {j = 1, ...,D) the components of the 
lattice position vector i in units of the lattice spacing. This charge configuration 
is exhibited in figure 4. For models with nearest-neighbor interaction, i.e. C/y = 
•^9+^ X]p '^i+pj with E denoting summation over nearest neighbors, the charge 
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Figure 4: ground state. 
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Figure 5: first excited state. 



configuration corresponding to the first excited state is given in fig. 5. Tlie energy 
of the charge distribution of fig. 5 is E[nl] — E[n^]+zd, where i?[n?], the ground 
state energy, is given by J2i i(l ~ ^^)- 

With the above values of E[n9] and E[nl] and keeping only the leading order 
term in Tc, eq.(43) becomes [see Appendix C] 



l-(l-z 




1- 



1 - (1 - z6 
1- (l + z6 



z 1- 



1 - (1 - z6 



1 - (1 - (z - 2)61)2 



(44) 

Reentrant behavior at low temperature occurs when the coefficient of the 
exponential is negative, namely when 



ai 



1- 



1- (l-z6 
1- (1 + ^6 



z l- 



l-{l-z6 



i-ii-iz-2)ey 



<0 



(45) 



In Appendix C we compute also the coefRcients 02 and 03 of the higher order 
exponentials in the expansion ([44|). In fig. 6 we plot the coefRcients oi, 02 and 
03 as a function of 9 for z = 6, i.e. for a 3-D array on a square lattice. One 
sees that the inequality ( [45| ) can be satisfied for values of 6 consistent with the 
physical constraint 6 < 1/z — 1/6. 

In fig. 7 we plot Tc versus a for 9 = 0.05 and z=6. In this plot we keep 
into account also the next two orders of (|4|) with coefficients 02 and 03. The 
resulting diagram exhibits reentrance in the insulating phase even for models 
with nearest neighbors interaction. 

In fig. 8 we plot ao = o:{Tc = 0) as a function of 9 for q integer and q half- 
integer and for z = 6. The plot shows that half-integer offset charges always 
favor superconductivity and that -at variance with the self-charging model- for 
non-diagonal interaction matrix there is always a range of a in which the system 
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Figure 6: Expansion cocfRcients of a as a function of nearest neighbor interac- 
tion 9. 
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Figure 7: Phase diagram for sniah critical temperatures with z = 6 and 9 = 0.05. 
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Figure 8: Broadening of the superconducting phase at T = with z = 6 and 
nearest neighbor interaction. 

is in the insulating phase. The plot also shows that for g/2e — 1/2 and T = 
the size of the superconducting region in the phase diagram depends on B. 

5 Discussion 

In this paper, using the path integral approach, we provided an explicit deriva- 
tion of the equation for the phase boundary for quantum Josephson junction 
arrays with offset charges and non-diagonal capacitance matrix. 

For the model with nearest neighbor capacitance matrix and uniform offset 
charge q — 1/2 (in units of 2e), using a procedure developed in Ig^li we were 
able to determine, in the low temperature expansion, the most relevant contri- 
butions to the equation for the phase boundary. For this purpose we explicitly 
constructed the charge distributions on the lattice corresponding to the lowest 
energies. 

Confirming the results of the numerical analysis of ref. S , we found a reen- 
trant behavior even with a short ranged interaction. Our analysis extends the 
results found in pfl to the situation in which offset charges are present and 
provides a physical picture of the states contributing to the reentrant behavior. 

For a model with diagonal capacitance matrix our analysis confirms the 
absence of reentrant behavior for the physical situation where the phase vari- 
able is 27r-periodic. The diagonal model with offset charge q = 1/2 exhibits 
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superconductivity for all the values of a = zEj/iEc, since in this case the su- 
perconducting order parameter is different from zero at zero temperature; this 
is evidenced by eq.(p3|) in Appendix A. An offset charge q = 1/2 tends to de- 
crease the charging energy and thus favors the superconducting behavior even 
for small Josephson energy Ej. 

A merit of the path integral approach, used in this paper, is that it allows 
to follow at each stage of the analysis the effects of the 27r-periodicity of the 
phase variable. In fact, one can untwist this periodicity by introducing a set of 
integers, so that the partition function factorizes as a product of a topological 
term, depending only on this set of integers, and a non-topological one explicitly 
evaluated in Appendix B. The Poisson resummation formula for the topological 
part of the partition function turns out very useful for the derivation of the low 
critical temperature expansion. 

It would be intersting to investigate the superconducting-insulating behavior 
in quantum JJA in lower dimensional models, where mean field theory is not 
expected to provide accurate results. For D = 1 there is evidence [0 for a new 
phase separating the superconducting and the insulating phase. The analysis 
of the phase diagram for this case should be carried out with different methods 
such as the renormalization group uK M . 

A Derivation of the self consistency equation 

m) 

With a uniform charge frustration q the pertinent Mathieu equation is given by 



dtf^ 2e dip \2e 



-^— 2 — 2i — -^ ^(~) — a < cos (/3 > cos (^ 



2 - p 

Mv)^^Mv) (46) 



Upon defining 

eq.(|46|) becomes 

^+(---cos^jp„ = (48) 

with A„ = En/Ec and v = —zEj < cos(p > /2Ec- Eq.(El) yields the canonical 
form of the Mathieu equations p5[ 



Mv) = e-'*^pn{v) (47) 



dx 



-j-^ + {X~'2vcos2x)y^0, (49) 



if one puts (p = 2x e i/jn — y- 
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The Mathieu equation has the well known periodic solutions [E5l 

even solutions with period tt 
with eigenvalues a2n{v) 



Ce2nix,v) 
Se2n+2{x,v) 



Ce2n+l{x,v) 
Se2n+l{x,v) 

71 = 0,1,2,.. 



odd solutions with period tt 
with eigenvalues b2n+2{v) 

even solutions with period 27r 
with eigenvalues a2n+i(w) 

odd solutions with period 27r 
with eigenvalues 62^+1 (w) 



If q/2e is integer, the periodic boundary conditions ipn{y^ = 0) = ipnif = 27r) 
singles out only the 27r-priodic Mathieu eigenfunctions ce2m se2n- With these 
eigenfunctions one may derive (M) lOJ. If q/2e is half- integer, the periodic 
boundary conditions together with (4S) single out the 7r-anti-periodic Mathieu 
eigenfunctions {i.e. pn is anti-periodic of 27r and periodic of Att). These are the 
Mathieu eigenfunctions ce2n+i and se2n+i- 

Since, near the critical temperature Tc, the order parameter < cosip > 
and V are small, apart from the phase factor e""''/^ (important only for the 
periodicity), to first order in v, eq.fiq) has the solutions 

^f = ;^(cos 



fcosf 



■ sm 






f) 



(2ri+l)y 



(2n-H)y 



(2Ti + 3)y 



4{2n+2) 

. (2,i + 3)y 



4(2n+2) 



cos 


(2 


n-1 


)v 




? 




"in 


8 

2ti 


n 
-1) 


»= -1 






2 





(50) 



8n 



(n = l,2,...) 
with the corresponding eigenvalues given by 

E! = Ecil + q) 



El = Ec{l-q) 



E: 



2n+l 



^2n+l 



Ec{2n+lf 



(51) 



(n = l,2,...) 

The expectation values of the superconducting order parameter on the eigen- 
functions ( pO[ ) are given by 

27r 
< V'nlcOS.^IV'n >= / d(fCOSip \ Ipnif) f ■ (52) 
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Using (bw, to the first order in v one gets 



<^?|cos(^|V?>=-i-f 

< V'l„+l|c0S(^|V'|„+i >= 8n(»+l) 

< V^2„+l|c0S(^|Vf„+i >= g;^^^ 
I (n = l,2,...) 



(53) 



Inserting (|5l|) and (|5^) in ([g) and keeping only the terms proportional to 
V ^< cos 93 >, one finds 



00 _i(2„+l)'i 

e y 



1 = a 



\y ^ 2>^ ^ 2n{n+l 



2e 



2Ee 

n=l 



-^(2«+l)^ 



(54) 



namely eq.(|ll| 



B The phase correlator 

In this appendix we want to elucidate the computation of the correlator defined 
in equation (Ea). For this purpose one should compute the path integral 



/i?0,e*W-^»(^')exp|/„''dr(-iCyi|)| 
/D</.iexp(/„^dr(-iajiii) 



(55) 



Fourier transforming (j)i{T) according to 

, +00 

n— — 00 

with < r < /3 and ujm — %""^i the numerator of d53) becomes 

/oo f ^ +00 

n ^'?^i^o n d4>i,nd(t)l^ exp i - -^ ^ ^ Cijw2 0i^„(/)?„ 

i n=l I- '^^ ij n=l 

n— 1 



(56) 



+ C.C. + - [(prfi - (psfi) 



n=l 

Upon integrating over the components c^r.Oj 0so one gets a factor (5r; 



11 / "<^i,0 



— 00 'J —00 



(57) 



(58) 
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where K is an irrelevant divergent constant which cancels against the denomi- 
nator. Using (|5^), ( p7| ) becomes 

OO n QQ I - 

• + 2^<Pi:no'5ri(e " -e " j-^^^LnOri^le " -6 " )l . 

The multiple Gaussian integral may be easily computed to give, up to an irrel- 
evant constant which cancels against the denominator, 

<5rs n n / ^<^i«'^'^in exp Y. 73'5ri(e^'^"^ - e^'^"^ )• 



i „=1 •'-'■- ^ ij 



e^C^r^ V^/1 — cosa;„(T — r'). 



'^ P n=l 



<5r.exp^~2e^C-/(|r-T'|-^^ ^'^' 



where — /3 < r — r' < /3. In the last step, the identity 



V.2 



CO 



x" v^^ 4 4cosijj„x, ^ ^ ,„„, 

has been used. This completes the proof of ( p6| ) 

C Low Tc expansion 

In this appendix we derive equation (H) and compute the next two orders whose 
coefficients are plotted in fig. 6. Using the notation (^1)' = (— 1)'^^ '"*", the 
ground state charge configuration n? can be written as 

The first excited states read 

(n^ + l)=nP(l-25ri); 

where the apex Ir means that this first excited state is obtained from the ground 
state by flipping the sign of the charge at the site r. Higher excitations may 
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be obtained from the ground state by flipping the sign of two charges at sites r 
and s and can be represented as 

The energy shifts are given by 

Ai - ^1 - s" = ^ ;7ir(-i)'-+'+i 



and 

A^" = £;[n?"] -E" ^ 2Ai + 2{-iy-^Urs . 

Note that, whereas the energy E^ of the charge configurations nl" does not 
depend on r, Eln^ ''"] depends on the relative position r — s of the charges whose 
sign has been flipped. 
Defining 



i?i^ = 



l-4Ej C/oj(nf + i)P 
1 



and 



i?^" = 



i-4Ej c/oj(nj- + i)F 
1 



l-4Ei Uoiir 



i)P' 



one may expand eq.(|4l|) for small critical temperatures (y oc Tc ^ 0), according 
to 



l+Er^' 



iA^ 



E 



iA^ 



ryis ' 



R" + Er R'^e-^^" + J2l^^ i?2,se-f A^.. ^ , 



1 



Eli '" . 4 A 1 



* P2r 



i?o 



E( 



i?^'ii'^ 



R' 



RO RO i?" 



)e « + 



(60) 



where Er^ts indicates a summation over pairs of different sites r, s, where each 
pair is counted only once. 

For a nearest neighbor interaction t/oj = Sqj + 9 ^ 6jp (where p denotes 
the vector connecting two neighboring sites) one has 

Ai = z9. 



A^" 



2(z-l) 

2ze 



r s = p 
r s 7^ p 
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R = r^—r 

1- (l-z6 






r^O,p 
r = 




r = p 


r 
r 


,s^O,p 
r^O,p 
= s = p 


r 


= p s = p' 



Substituting these relations in (pOD, one obtains the expansion for sniaU tem- 
peratures of the critical line equation, up to the first four orders 

a=fl-il- ze)A xfl + a, e-^''" + a^ e^^^'-'^' + a^ e~iA . (61) 

oi is given in (pa), 02 is equal to 



^ '[ l-il-iz-2)9yj^\' l-(l + (z- 2)0)2 

and as is given by 

l-{l-z0fY fl-{l-z9)^\ ( l-{\-zO^^ 



\-(i + zeY) \i-(i + zef)^^^^ ^^ \i - {\ - {z - 2)eY ^ 

2 1 - (1 - zef ( 1 - (1 - z9f \ 1 - (1 - zef ( 1 - (1 - z9f 



"^^ 'l-(l- (z- 2)61)2 (^l-(l-(z -2)6')2 ^)^'''i-{i + zeY\l-{l-{z-2)eY ^ 

,^. i-{i-zef f i-{i-zef ^V ^(^-i) A 1 - (1 - zef 



l-{l-{z-2)6f) \l-{l + zeY J 2 \ l-(l-(z- 4)0)2 

The condition for the reentrant behavior is ai < 0. In fig. 6 we plot the coefh- 
cients ai, 02, as as a function of 9. In fig. 7 we plot the critical equation (pi]) 
with e = 0.05 and z = &. 
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